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Abstract 

The bosonization and duality rules in three-dimensions are applied to analyze some fea- 
tures of superfluids and superconductors. The energy of an ensemble of vortices in a su- 
perfluid is recovered by means of a kind of bound which, to some extent, shares similarity 
with the Bogomol'nyi bound. In the case of superconductors, after recasting the partition 
function in the form of a pure effective gauge theory, the existence of finite energy vortex 
solutions is discussed. 
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1 Introduction 



In the last years many efforts have been undertaken to generaUze the bosonization techniques 
^ to fermionic models in (2 + 1) dimensions |1]-|HI- Although a complete framework is still 
lacking, these efforts have led to a set of interesting and promising results, successfully applied 
to describe transport properties of interacting electronic systems, such as the universal behavior 
of the Hall conductance The bosonization rules so far obtained are efficiently derived by 
means of the path integral 0-13) being summarized as follows. 

Let jfj_ = ipj^ip be the U{1) conserved current of a fermionic system in (2 + l)d. It turns out 
that jfj, is mapped into a topologically conserved current [2]-[H], according to 

= V'7m'0 — ' = (^/luxduAx . (1.1) 

From this equation we see that the bosonization of the fermionic current is achieved through 
the introduction of an Abelian gauge field Ax. Eg. 1)1.1(1 generalizes the well known bosonization 
rule in (1 + l)d, as expressed by the equivalence between the massive Thirring model and the 
sine-Gordon theory The meaning of the equation (|1.1|) can be stated in a more precise way 
at the level of the current correlation functions, namely 

{jlJ.l---jflJSF = {j'lli---jJ^n)sB[A] ) (1-2) 

where Sp is the free massive Dirac action 

= / c?3xV5(i7''^ + m)V , (1.3) 

and 5'b[^] is the corresponding bosonized action, which displays the important property of 
being gauge invariant. We underline that, in general, a closed form for S^i^] is not available, 
since its evaluation requires the exact knowledge of the fermionic determinant in (2 + l)d for 
an arbitrary gauge field configuration. However, relying only on gauge invariance, one can show 
[3 IH] that 5b [^] consists of the Chern-Simons action with the addition of an infinite series of 
terms, depending on the curvature = ^^/luxduAx, i.e. 

Sb[A] = ^^Scs[A] + R[F], (1.4) 
47r \m\ 

where 5'c5[^] is the Chern-Simons action 

Scs[A] = ^Jd^xe^^xA^d^Ax, (1.5) 

and R[F] denotes the remainder, higher order contributions in the curvature F^ |21-[ZI- The 
presence of the Chern-Simons term in the expression ()1.4() is expected |11| since the fermionic 
action ()1.3|) contains the parity breaking mass term m^/^^/^ . 

It is worth mentioning that the bosonization rule (jl.2() has been proven to hold also in 
the presence of a large class of interaction terms I[j^] depending only on the fermionic current, 
namely 



SF[i^,ij]+I[f] ^ SBlAl+Iie^'^Pd.Ap], (1.6) 
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Notice that the bosonization of the interaction term is obtained by replacing the fermionic 
current j'^ = V'T^tV'by the topological current = e^y\dyA\, without changing the functional 
form of the interaction term This important property gives a universal character to 

the bosonization rules H1.6|) . For instance, in the massive Thirring model with a non- 

polynomial current-current interaction has been bosonized into the (2 -|- l)d Born-Infeld gauge 
theory, according to eas. ()1.6() . 

More recently, the authors mj have been able to generalize this construction, originally used 
for the bosonization of fermionic actions, to the case of an arbitrary field theory model with a 
conserved global U (1) charge. It has been proven that, regardless of the spin and statistics of the 
fields in the original action, any three-dimensional model with an Abelian conserved charge can 
be mapped into a dual Abelian gauge theory jl4j . This mapping acquires therefore the meaning 
of a duality, in view of its more general defining characteristics. As the most distinctive feature 
of this mapping, we mention that the conserved Noether current j^^^^ corresponding to the U{1) 
global symmetry has, again, a topological current as its dual partner in the Abelian gauge 
theory, i.e. 

i^^'^ jj = e^uxd.Ay^ . (1.7) 
Equations ()1.6|1 have also their generalization, which reads 

S[^, (/>t] + ^ Sduail^] + Ae^^uxdM , (1.8) 

where 4>,4>^ are a pair of generic complex U{1) charged fields, S[(p,(l)^ is the corresponding 
invariant action, and collects the interaction terms depending only on the Noether 

current. 

In the present work we shall apply the dual mapping (|1.7|) — (|1.8j) to analyze some features 
of superfluids and superconductors, the aim being to show that the dual description in terms of 
gauge theories can be effectively worked out. As first example, we shall consider the superfluid 
'^He system in the presence of an ensemble of vortices. The dual mapping will be used to evaluate 
the energy of the vortices. In particular, we shall be able to prove that the energy configuration 
can be obtained by means of a relationship which, to some extent, shares close similarity with 
the Bogomol'nyi bound of superconductors In this case the mapping (|1.7|) — (|1.8|) will have 
the meaning of a duality between vortices and electric charges. A vortex line can be regarded 
in fact as the world-line of a particle whose electric charge is proportional to the topological 
integer corresponding to the so called strength of the vortex. Although the duality between 
vortices in ^He and electric charges is already known |151 116j . we point out that it is a natural 
consequence of the relation (|1.7|) . In the second example, after recasting the partition function 
of the Landau-Ginzburg model for a superconductor in the form a pure three-dimensional gauge 
theory, we shall analyze the existence of vortex configurations in the resulting effective theory. 
We remark here that these examples will be worked out in the Euclidean three-dimensional 
space rather than in (2 -|- l)d. In this case, the bosonization rules (|1.6|) and the dual mapping 
(|1.7|) — H1.8() turn out 14 to describe the static properties of three-dimensional systems. 

It is worth underlining that the possibility of describing three-dimensional U{1) charged 
models in terms of the corresponding dual gauge theory can be useful in order to improve our 
understanding of the nonperturbative aspects of these models. For instance, in the case of 
superconductors, it turns out that the resulting effective dual gauge theory is described by an 
action containing a term of the type / d^x {Fij -\- Jij)'^, where Jij accounts for the presence of 
vortex lines. This kind of action is deeply related to the formulation of the so called compact 
U{1) gauge theories in continuum space-time. These models are object of intensive investigation 
in order to obtain a better understanding of the phase structure of confining theories and in 
order to test new order parameters, as in the case of the operator ^^i^ recently proposed in the 
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context of three-dimensional compact QED Also, this kind of action naturally arises in 

SU{2) Yang-Mills theories which, after performing 't Hooft abelian projection, turn out to be 
described by an effective U{1) compact theory which is expected to be the basic tool for 
the dual superconductivity mechanism for the confinement of quarks and gluons. 



2 Vortices in the ^He Superfiuid 

2.1 Generalities 

In this section we shall briefly review some basic features of the ^He superfiuid in the presence 
of an ensemble of vortices |15j . Let us begin with the expression of the Landau-Ginzburg free 
energy for the ^He in the superfiuid phase 

E = [ d^xdiip+di^ + ^ -a^Y , (2.9) 
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where the complex field (p stands for the superfiuid order parameter and the index i runs from 1 
to 3. For temperatures T <^T\ well below the critical value Tx — 2.\^K corresponding to the so 
called superfiuid A— transition, only the phase fluctuations of the order parameter are relevant, 
namely 

(^(x) = f7e^^(^) , (2.10) 

where the parameter a is related to the density of the superfiuid component ^Hl) identified with 
the Bose-Einstein condensate. Therefore, for the Landau-Ginzburg free energy we get 



E = a^ J d^xdada . (2.11) 



The phase configurations which minimize the functional ()2.11() are those obeying the Poisson 
equation 

dida = . (2.12) 

The presence of vortices in the superfiuid is associated with the existence of lines of singularity in 
the phase 7, which turns out to be multi-valued. These lines can be of any shape and they never 
end inside the superfiuid. They can be closed or run to infinity^. For instance, an ensemble of 
N vortices corresponding to a set of closed lines I = ^1 •••j-^j is described by 

eijkdjdk-f{x) = 2vrn(^) J^^ dyi5\x - y) , (2.13) 

where denotes a smooth closed curve and n^^^ is an integer. From this expression it is 
manifest that the phase 7 is singular along each curve , the singularity being expressed by the 
line integral of the 5-function, i.e. J^i dyi5^{x — y). Since the curves |-^^| are closed, it follows 
that 

eijk^^^J^kl{x) = . (2.14) 



^In this last case they extend from —00 to +00, being meant to be closed at infinity. 
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Also, from eq. (j2.1Hj) . we see that to each vortex Une there corresponds an integer rS-^\ known 
as the strength of the vortex. Before ending this short summary, it is useful to mention that the 
vortices have a finite size, looking alike thin tubes of radius where ^ is the coherence length. 
The region r < ^ is called the core of the vortex. The superfluid does not penetrate inside the 
core region, which contributes to total energy of the fluid. For an estimate of the core energy 
see [T^ . 



3 Evaluation of the energy of an ensemble of vortices 

Let us turn now to the evaluation of the energy E in the expression ()2.11() for the ensemble of 
vortices specified by the eq. (|2.13|) . In order to make use of the dual mapping (|1.7j) — (|1.8|) , we 
begin by computing the U{1) current j^^^^ for the Landau-Ginzburg functional (|2.9j) . We easily 
find 



jf = I V - V^) = <^'da . (3.15) 

According to eqs. (|1.7|) . (jl.8|) . the free energy (j2.11|) has a dual description in terms of a gauge 
theory, the dual gauge field Ai being identified through the topological current jf = eijkdjAk, 
namely 

a'^da = CijkdjAk . (3.16) 
As it will become clear in the following, it is useful to introduce the quantity Q defined by 




(3.17) 



where, for the time being, Ai denotes a generic field configuration, not yet identified with the 
dual gauge field of the eq. ()3.16() . Obviously, we have 



e > , (3.18) 

the equality holding when Ai is the dual field. The condition ()3.18|) can be seen as a kind of 
bound which, to some extent, shares similarity with the Bogomol'nyi bound of superconductors 
jl7j . The bound is saturated when G = 0, i.e. when the gauge field is identified with the dual 
field of eq . (13.161) . Notice also that, as in the case of the Bogomol'nyi bound, the expression 
(|3.16|) is of first order in the derivatives of the fields. 

Expression (|3.17|) is easily worked out, yielding 

e = Jd^x (^a^a,75a + l^.^Fij - 2ahijkAdjdk7^ . (3.19) 

Requiring now that Ai is precisely the dual gauge field of eq. ()3.16|) , the bound gets saturated so 
that G = 0. Thus, for the energy of the ensemble of vortices we get 



E = 'J^ J d^x^a^^J = J d^x {^-^PijF^i + 2^i^i) > 
with Fij = {diAj — djAi) and 

/,(x) =5]27rn(^) / dy,6^{x-y) . 



(3.20) 



(3.21) 
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The equation (|H.2()|) shows that, according to the dual mapping p.7j) — (|1.8|) . the energy of 
an ensemble of vortex lines can be obtained by evaluating the action of an electromagnetic 
field interacting with the external current Ii{x). Furthermore, we notice that the current Ii{x) 
in ea. ()3.21|) corresponds to a system of charged particles whose electric charges are given by 
= 4Trn^^\ In addition, the vortex line integrals dyi6^{x — y) can be associated to the 
world-lines of the particles, making manifest the duality aspects between vortices in the ^He 
superfluid and a system of electric charges. We also remark that, from the Poisson equation 
(Einj), it follows that 

eijkdidjAk = , (3.22) 

which implies the absence of magnetic charges. 

In order to evaluate the electromagnetic action in ea. ()3.2(J| ). we act with the operator emnidn 
on both sides of the equation (|3.16j) . obtaining 

dkdkA, = -a^Ii , (3.23) 

where use has been made of the gauge condition dkAk = 0. For the dual gauge field we get 

Ai{x) = ^ / d'y—^ L{y) . (3.24) 
47r J \x — y \ 

Thus, 

Finally, the energy of the ensemble of vortices if found to be 

E = f d^xdndn = a\y n^^'^n^-^'^ f dlt^ ■ f d^^ ^- , (3.26) 



i,j 



which coincides in fact with the well known expression given in |151 116j . We have thus recovered 
the energy of an ensembles of vortex lines by using the dual mapping. 



4 Effective gauge theory for superconductors 

In this second example we shall consider the partition function of a superconductor. After 
recasting it in the form of a pure gauge theory, we shall discuss the existence of finite energy 
vortex solutions in the resulting effective theory. 

The partition function of the Landau-Ginzburg free energy of a superconductor reads 

Z=l DAiDip+Dipexp - J d^x (^^F.^Fij + {Vi^)+ {V,^) + V{^) ^ , (4.27) 

where ip stands now for the order parameter describing the condensation of the Cooper pairs 
and 

V{ip)=^[^+^-a^)' . (4.28) 
The covariant derivative is defined as 

Vi = di- iqAi 
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with q = 2e the charge of a Cooper pair. In order to obtain a representation of the partition 
function H4.27() in terms of an effective gauge theory, we first write the order parameter ip as 

^ = {a + p)e^^(") , (4.29) 

where p accounts for the fluctuations above the condensate a. Further, we decompose the phase 
9 according to 

9 = 9' + 9\ (4.30) 
where 9', 9^ denote the regular and the singular part of 9 '20', namely 



eijkdjdk9\x) = , 

eijkdjdk9\x) = I^{x)=2Tm f dy^d^{x-y) , (4.31) 

where n is an integer. As in the case of the ^He superfluid, the quantity Ii{x) specifies the 
singularity of 9^{x), which turns out to correspond to Abrikosov-Nielsen-Olesen (ANO) [2J| flux 
tubes inside the superconductor. We shall limit here to consider the case of a unique smooth 
curve L, the example being easily generalized to an ensemble of curves. 

Making use of the change of variables 



A, ^ A, + -{di9' + di9') , (4.32) 



and noticing that, due to ea. (|4.31() . 



for the partition function we get 

Z=N j DA^D9'e-^-«^^'-^'> , (4.33) 

where is a normalization factor and £^efr(^)<^) stands for the effective energy functional 
corresponding to 

E,s{A, J) = j (fx Q {Fi, + J,,f + a^q^A^^ + E{A) , (4.34) 

E[A) denoting the contribution coming from the integration over the radial fluctuations p, 
namely 



e 



EiA) ^ I Dpexp - J Sx (dipdip + q^p{p + 2a) A^ + V{p) 



(4.35) 



Although we shall be interested in the study of the classical vortex solutions of £'efi(^; J)-, it is 
worth mentioning that the integration over the singular variable D9^ in the expression ()4.33|1 
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can be specified in a more precise way as an integration over all possible smooth curves {L}, as 
discussed in ^1 120] . 

The expression (j4.33() yields the representation of the partition function Z in terms of the ef- 
fective gauge theory we were looking for. The resulting effective functional Ecq^A, J) in eg. (|4 .3411 
is sometimes referred as to a compact gauge theory, due to the presence of singular configura- 
tions, encoded in the line integral Jij and in the functional measure D^'^ llSl 120] . 

Let us turn now to the study of the existence of finite energy vortex configurations for the 
resulting energy functional Ecsi^, J)- We shall begin by working out the so called London limit, 
corresponding to A — > oo. 

4.1 Vortex solution in the London limit 

In the London limit, i.e. A ^ oo, the radial part p of the order parameter gets frozen 

<f = cje*^(^') . (4.36) 
As a consequence E{A) = 0, so that Ees{A, J) reduces to 

^Lond(A, J) = J d^x Q {Fi, + J,,f + ^m^AA,^ , (4.37) 

with 

= 2a^q^ . (4.38) 
The configurations which minimize the functional -ELond(^) J) those for which 



diFij + m^Aj - diJij = . (4.39) 



From ea. ()4.39() it follows that Ai is transverse 

diAi = . (4.40) 



Therefore, expression (|4.39|) becomes 



From 



{d^ - m^)A, = ^ei.kdjh • (4.41) 



{-d'^ + m^)g{x -y) = 6^{x-y), (4.42) 

(— a^+m^j 4tt \x — y\ 



we get 



Ai{x) = -^KjkdJ I tfy-^^e-^^^^-y^Uy) , (4.43) 
Att q ■'J F — 2/1 

which, from ea. (|4.31|) . takes the final form 
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Mx) = -^ei.kd] I dyfc^^e-'"!--^! . (4.44) 

This expression yields the solution of the equations of motion for a generic smooth curve L. It 
corresponds in fact to a vortex solution localized around the curve L, on a scale of the order 
1/m. For a better understanding of expression H4.44() let us work out in details the case in which 
L is a straight line coinciding with the z-axis. In this case, the line integral can be evaluated in 
closed form, giving 



Tl 

Ai = eijsdjKoimx^) , (4.45) 

q 

Tl Tl 

Ax = --dyKo{mx±) , Ay = -dxKo{mx±) , A^ = , 



where Kq is the Bessel function and x± = ^/aP + y^ is the distance from the origin in the 
(x, y) —plane. The magnetic field B = is found 

71 1 

B = d^Ay - dyA, = {x±d.,^Ko{mx±)) . (4.46) 

Also, for the flux of B through the circle at infinity in the (x, y) —plane, one obtains 

$ = / (fxB = / dx^d^^ {x±d.^^Koimx±)) = . (4.47) 

Jsl, q Jo q 

The gauge configuration H4.45() is nothing but the London limit of the (ANO) vortex solution. 
In particular, the equation ()4.47() expresses the well known flux quantization of the ANO vortex 
in superconductors. We have thus checked that the effective gauge functional £^Lond(^)<^) of 
eq. (|4.H7)) admits a vortex solution with the correct quantization of the magnetic flux. Moreover, 
because of the limit A ^ oo, the configuration H4.45() does not fulfill the requirement of finite 
energy. As we shall discuss in the next section, this feature can be handled by going beyond the 
London limit. 



4.2 Finite energy solutions 

So far we have seen that the functional Ei^ond{A, J) of ea. ()4.37|) admits vortex configurations, 
although they do not meet the requirement of finite energy. This means that the classical 
configuration ()4.45|1 gives a divergent value for £^Lond(^) J)- Indeed, from ea. (|4.37|l we have 

SLond(^, J)= J d^x Q^i (-a^ + m^) Ai + ^A,eijkdjh + ^44^ ■ (4.48) 



A = -z . .7^^2, 9jh, (4.49) 



From 

q (—9^ + m 
one easily gets 

E^ondiA, J) = ^j ^'^4p^2^^4 . (4.50) 
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In the case in which Ik is a straight hne coinciding with the z-axis, expression (|4.5())) is easily 
written in momentum space, yielding 



""^"^^ - ^ d\^-^ , (4.51) 



L Iq^ J k\ + 



m 



where Ei^o^^/L is the energy per unit length^ and k ± is the two-dimensional vector k± = 
{kx,ky,0). We see that E'Lond/-^^ displays an ultraviolet logarithmic divergence in momentum 
space. We need therefore to go beyond the London limit, which amounts to evaluate the quantity 
E{A) in expression ()4.34() . This is a rather difficult task, since it requires the evaluation of the 
functional integral of ea. (|4.35|) . whose exact form is unknown, even if the cubic and quartic terms 
contained in the potential V{p) are neglected^. Nevertheless, due to the gauge invariance, we 
expect that E{A) will contain nonlocal terms in the curvature Fy. It is worth to point out here 
that these nonlocal terms naturally arise in the evaluation of three-dimensional effective action 
as, for instance, in the case of the bosonized action corresponding to the fermionic determinant 
[7||n]. Moreover, as shown in these nonlocal terms provide a useful ultraviolet regularization, 
allowing in fact to fulfill the requirement of finite energy. Following [22] , a simple understanding 
of the regularizing mechanism due to the presence of nonlocal terms can be achieved by replacing 
in the effective functional ii^Lond(^5 J) of eq. (|4.48j) the field strength Fij by 

Fij ^ d{d)Fi,j , (4.52) 

where O is a suitable nonlocal operator associated with a kernel 0{x — y), according to 

O/] (x) = J d%0{x - y)f{y) . (4.53) 

With the replacement (|4.52l) . the functional -ELond(^) J) becomes 

i^Lond(AJ) ^i?0(AJ) , (4.54) 



with 



Eo{A, J) = j (^^FijO^Fij + ^FijOJij + ^JtjJ^j + ^m^AiA,^ (4.55) 
= J d^x Qyli (-O^a^ + m^) Ai + ^AiCijkOdjlk + ^hh 



For the vortex solution we get now 



from which it follows that 



where 



= — I n2%\ 2, O0jh , (4.56) 



fA, ,,'^ (4,57) 



L 2g2 J ^klO^{k^) + m^ ' 



0{p) = I d^xe-'P^'Oix) , (4.58) 



^As usual, we have factorized the integration over the 2:-axis, i.e. J dz = L, due to the independence of the 
integrand in eg. 14. 50(1 from the z coordinate. 

^See \2'2\ for a recent attempt to go beyond the London limit in the case of the four-dimensional euclidean dual 
Abelian Higgs model. 
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is the Fourier transform of the kernel 0{x). In the local case, i.e. 0{k±) = 1, we recover the 
expression ()4.51() with its ultraviolet logarithmic divergence. However, in the case where in the 
ultraviolet limit 0{k±) behaves as fc^, a > 0, the energy per unit length is finite, no matter 
how small a is. We see thus that the requirement of finite energy can be fulfilled by means of 
the inclusion of suitable nonlocal gauge invariant terms. These terms are expected to show up 
as the result of the functional integration over the radial fluctuations of the order parameter in 
en. lfOnl) . 

5 Conclusion 

In this paper the bosonization rules (|1.6|) and the dual mapping (|1.7|) — (|1.8j) have been applied 
to analyze some features of superfluids and superconductors. In the case of superfluids, we 
have been able to recover the energy of an ensemble of vortex lines in the ^He superfluid by 
use of the dual mapping (|1.7|) — (|1.8|) . In particular, we have derived a kind of bound for the 
energy configuration which, to some extent, shares similarity with the Bogomol'nyi bound of 
superconductors. When the bound is saturated, the dual mapping applies and one recovers 
the energy of an ensemble of vortex lines in terms of a dual electromagnetic action for charged 
point-like particles. The world-lines of these charges correspond to the vortex lines. In the case 
of superconductors, after recasting the partition function in terms of an effective pure gauge 
theory, the existence of finite energy vortex solutions with the correct flux quantization has 
been established. 
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